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ABSTRACT 

Measuring the intrinsic shape and orientation of dark matter (DM) and intracluster 
(IC) gas in galaxy clusters is crucial to constraining their formation and evolution, 
and for enhancing the use of clusters as more precise cosmological probes. Extending 
| our previous works, we present for the first time results from a triaxial joint analysis 

of the galaxy cluster Abell 1835, by means of X-ray, strong lensing (SL) and Sunyaev 
Zel'dovich (SZ) data. We parametrically reconstruct the full three-dimensional struc- 
ture (triaxial shape and principal axis orientation) of both the DM and the IC gas, and 
psj ■ the level of non-thermal pressure of the IC gas. We find that the intermediate-major 

1 and minor-major axis ratios of the DM are 0.71±0.08 and 0.59±0.05, respectively, and 

ON . the major axis of the DM halo is inclined with respect to the line of sight at 18.3 ± 5.2 

deg. We present the first observational measurement of the non-thermal pressure out 
to -R200 * which has been evaluated to be a few percent of the total energy budget in the 
. internal regions, while reaching approximately 20% in the outer volumes. We discuss 

the implications of our method for the viability of the CDM scenario, focusing on the 
concentration parameter C and the inner slope of the DM 7 in order to test the cold 
dark matter (CDM) paradigm for structure formation: we measure 7 = 1.01 ± 0.06 
and C = 4.32 ± 0.44, values which are close to the predictions of the CDM model. 
The combination of X-ray/SL data at high spatial resolution, capable of resolving the 
cluster core, with the SZ data, which are more sensitive to the cluster outer volume, 
' allows us to characterize the level and the gradient of the gas entropy distribution 

J-j . and non-thermal pressure out to i?200, breaking the degeneracy among the physical 

models describing the thermal history of the ICM. 

Key words: cosmology: observations - galaxies: clusters: general - galaxies: clusters: 
individual (Abell 1835) - gravitational lensing: strong - X-rays: galaxies: clusters - 
cosmic microwave background 



1 INTRODUCTION 

The cold dark matter (CDM) paradigm has been remarkably 
successful at predicting the large-scale distribution of mat- 
ter in the Universe as well as its observed evolution from the 
earliest epochs to the present day. A fundamental prediction 
of N-body simulations is that the CDM halos follow a self- 
similar density profile, with the logarithmic slope of the DM 
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7 following a shallow power l aw at small radii (7 ~ 1), which 
then steepens at larger radii jNavarro et al.|[l996l ). However, 
a comprehensive physical explanation for the origin of such 
a profile is still lacking. Moreover, in recent years much in- 
terest has been shown in possible discrepancies that remain 
for th e observed and predicted inner density profile of struc- 
tures (|Sand et al.ll2008l ; lLimousin et al.ll2008l : lNewman et al.1 

In this perspective, clusters are an optimal place to 
test the predictions of cosmological simulations regarding 
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the mass profile of dark halos, and to in general cast light 
on the viability of the standard cosmological framework 
consisting of a cosmological constant and cold dark mat- 
ter (ACDM) with Gaussian initial conditions, by comparing 
the measured and the predicted physical parameters (e.g. 
concentration parameter, inner slope of the DM). For exam- 
ple, observations based on the combination of strong lens- 
ing and stellar kinematics yielded flat inner slopes 7 ~ 0.5 
for two well-stu died clusters (MS21 3 7-23 and Abell 383, 
ISand et al.ll2008h . and lNewman et all (l201ll ) derived a shal- 
low cusp with 7 < 0.3 (68%) for Abell 611, raising doubts 
on the predictions of the CDM scenario. Other studies lead 
to large scatter in the value of 7 from one cluster to an- 
other, but these determinations customarily rely on the stan- 
dard spherical modeling of galaxy clusters. Possible elonga- 
tion/flattening of the sources along the line of sight, as well 
as the degeneracy of 7 with other parameters, i.e. the con- 
centration parameter and the scale radius, likely affect the 
estimated values of 7 l|Morandi et aljfeoiol ). 

Clusters are also an optimal tool to constrain the cos- 
mological parameters provided that we can accurately de- 
termine their mass. For example, comparison of the cluster 
baryon fractions /(, to the cosmic baryon fraction can pro- 
vide a direct constraint on the mean mass density of the 
Universe, Q m (|Ettori et alj|2009h . while the evolution of the 
cluster mass function can tightly constr ain Oa and the dar k 
energy equation of state parameter w (M antz et alj 12010). 
Cluster mass profiles can be probed through several indepen- 
dent techniques, relying on different physical mechanisms 
and requiring different assumptions. 

So far, analysis of the cluster X-ray/SZ emission and of 
the gravitational lensing effect are among the most promis- 
ing techniques to estimate galaxy cluster masses. Concern- 
ing the former, the cluster mass can be measured by study- 
ing the IC gas emission under the ass umption of hydro- 
static equilibrium (HE, see lSarazi n Hl988T ). Indeed the IC gas 
emits via both thermal bremsstrahlung in the X-ray band 
and inverse Compton with the photons of the CMB spec- 
trum , a process known as the S unyaev-Zel'dovich (SZ) ef- 
fect (|Sunvaev fc Zeldovich|[l970l) . The X-ray mass estimate 
is less biased compared with SZ and lensing-derived masses 
with regard to projection effects, because the emission is 
traced by the square of the gas density. 

The advantage of the SZ effect compared to the X-ray 
emission is the possibility of exploring clusters at higher red- 
shift, because of the absence of the cosmological dimming. 
Moreover, since the SZ intensity depends linearly on the den- 
sity, unlike the density squared dependence of the X-ray flux, 
with the SZ effect it is possible to study clusters without the 
systematic errors caused by the presence of sub-clumps and 
gas in multi-phase state and to study the physics of the ICM 
well beyond the regions constrained with X-ray observations 
« 0.3 -0.5 #200). 

On the other hand, the gravitational lensing effect al- 
lows for the determination of the projected surface mass den- 
sity of the lens, regardless of its dynamical state and inde- 
pend ent of the assumption of HE ( Miralda- Escude fc Babuj 
Il995l ). Unfortunately, in most cosmological applications the 
projected mass is not the interesting quantity. Rather, we 
need to measure the three-dimensional mass profile, custom- 
arily by assuming spherical symmetry. Lensing mass mea- 



surements are also appreciably prone to contamination from 
foreground and background sources. 

Knowledge of the intrinsic shape and orientation of ha- 
los is crucial in order to obtain unbiased determinations 
of their masses, inner slope of the DM and concentration 
parameter via e.g. X-ray, SZ and lensing data. From this 
perspective, clusters are commonly modeled as spherical 
systems whose intracluster (IC) gas is in strict HE (i.e., 
the equilibrium gas pressure is provided entirely by ther- 
mal pressure), assumptions that are only rough approxima- 
tions, leading to large biases in the determination of the 
cluster mass and hence on the desired cosmological param- 
eters. Indeed N-body simulations indicate that DM halos 
are triaxial with intermediate-major and minor-intermediate 
axis ratios typically o f the order of ~ 0.8 (|Shaw et al.ll2006l ; 
IWang fc Whitdl2009i ) , while hydrodynamical numerical sim- 
ulations suggest that the plasma in apparently relaxed sys- 
tems may also be affected by additional non-equilibrium 
processes, which serve to boost the total pressure and there- 
fore cause an undere stimate of the clust e r mass from X- 
rav/SZ observations (lAmeglio et ail 120091 ; lLau et all 120091 ; 
iMeneghetti et alj|2010h . 

On the observational side, only a few works have 
tried to infer shape or orientation of single objects 
dOguri et all 120051; ICorless et all 120091 ; iMahdavi fc Chand 



l201ll ; IMorandi fc LimousinllOlld), and the non-thermal pres- 
sure support llMahdavi et all 120081; ISanders et al.l l20ld : 



I Richard et a.1.1 [20101 : Morandi fc Limousin Olid). By means 
of a j oin t X-ray and lensing analysis , iMorandi et al.l (|2010l . 
1 11 al lbh : IMorandi fc Limousin! (|01ld ) overcame the limita- 
tion of the standard spherical modelling and strict HE as- 
sumption, in order to infer the desired three-dimensional 
shape and physical properties of galaxy clusters in a bias- 
free way. A triaxial joint analysis relying on independent 
and multifrequency datasets for galaxy clusters can relax 
the assumptions customarily adopted in the cluster analysis 
and give us additional insights into the underlying physics 
of these objects. 

Extending the findings of our previous works, in the 
present paper we recover the full triaxiality of both the DM 
and the ICM, i.e. ellipsoidal shape and principal axis orienta- 
tion, and the level and the gradient of non-thermal pressure 
for the galaxy cluster Abell 1835. This cluster is a lumi- 
nous cool-core galaxy cluster at z = 0.253 and it is an opti- 
mal candidate for a triaxial joint analysis via X-ray, SZ and 
lensing techniques, because of its very relaxed dynamical 
appearance and its exceptional strong lensing system. We 
discuss the implications of our findings for the viability of 
the CDM scenario, focusing on the concentration parameter 
and inner slope of the DM. 

The availability of SZ data out to .R200 allows us to 
infer the properties of the ICM in the outskirts of the galaxy 
cluster. An accurate measurement of the properties of galaxy 
clusters out to large radii provides critical insights into the 
physics of the ICM and offers a direct probe of the assembly 
history of structure formation on the largest scales; it also 
enhances the use of clusters as cosmological probes, since the 
physics of the IC gas in the outer volumes is relatively simple 
and nearly self-similar. In particular we compare our findings 
with the results of hydrodynamical numerical simulations for 
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the density, temperature, entropy and non-thermal pressure 
profiles out to the virial radiu^J. 

Throughout this work we will assume a flat ACDM 
cosmology: the matter density parameter will have the 
value f2 m =0.3, the cosmological constant density param- 
eter £2a=0.7, and the Hubble constant will be Ho = 
70 km s _1 Mpc _1 . At the cluster redshift and for the as- 
sumed cosmological parameters, 1 arcsec is equivalent to 3.9 
kpc. Unless otherwise stated, quoted errors are at the 68.3% 
confidence level. 



2 STRONG LENSING MODELING 

A detailed mass model of Abell 1835 and a description of the 
data used will be presented in a forthcoming publication. In 
this Section, we summarize it. 



2.1 Multiple images 

T he last strong l e nsing model of Abell 1835 was presented 
bv lRichard et all |201oh . 

Since then, Abell 1835 was observed with WFC 3 on- 
board the Hubble Space Telescope (HST) in the F110W 
and F160W filters (Program 10591, PI: Kneib). These new 
imaging data allowed the identification of new multiply im- 
aged systems in the cluster core. Moreover, a spectroscopic 
campaign targeting multiple images was carried out using 
FORS 2 on the VLT (Program 087. A-0326, PI: Julio) and 
yielded a spectroscopic measurement for some multiple im- 
ages. 

The model presented bv lRichard et al l (|2010h was based 
on 7 multiply imaged systems, 2 of them having a spectro- 
scopic redshifts measured. In this paper, the model is based 
on 8 multiply imaged systems, 6 of them being spectroscop- 
ically confirmed. For the remaining systems, the redshifts 
will be let free during the optimization. These images are 
reported in Table [1] and shown on Fig. [1] 



2.2 Mass Distribution 

The model of the cluster mass distribution comprises three 
mass components described using a d ual Pseudo Isothermal 
Elliptical Mass Distri bution (dPIE, iLimousin et al.1 120051 : 
lEKasdottir et al1l2007l ). parameterized by a fiducial veloc- 
ity dispersion a, a core radius r corc and a scale radius r s : 
(i) a cluster scale dark matter halo; (ii) the stellar mass in 
the BCG; (iii) the cluster galaxies re presenting local pertur - 
bation. As in earlier works (see, e. o. ILimousin et al.1 007a1. 
empirical relations (without any scatter) are used to relate 
their dynamical dPIE parameters (central velocity disper- 
sion and scale radius) to their luminosity (the core radius 
being set to a vanishing value, 0.05 kpc) , whereas all geomet- 
rical parameters (centre, ellipticity and position angle) are 
set to the values measured from t he light distribution. B eing 
close to multiple images, following [Richard et al.1 ( |2010h . two 
cluster galaxies are modeled individually, namely PI and P2 
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Table 1. Multiply imaged systems considered in this work. Red- 
shift measurement and estimation will be presented in a forth- 
coming publication. 



(Fig[T}. Their scale radius and velocity dispersion are opti- 
mized individually. We allow the velocity dispersion of clus- 
ter galaxies to vary between 100 and 250 kms , whereas 
the scale radius was forced to be less than 70 kpc in order 
to accou nt for tidal stripping of their dark matter haloes 
Limousin et"alTl007bi . 120091 : iNataraian et al.ll2009l : 



1 Hereafter we equate the virial radius with i?200i the radius 
within which the mean total density is 200 times the critical den- 
sity of the Universe at the redshift of the cluster. 



IWetzel fc Whitdl201Cl . and references therein.) 

Concerning the cluster scale dark matter halo, we set 
its scale radius to 1 000 kpc since we do not have data to 
constrain this parameter. 

The optimization is p erformed in the im age plane, using 
the LenstooiQ software \ Julio et al.ll2007f ). 



2.3 SL results 

The RMS in the image plane is equal to 1.4 . In good agree- 
ment with lRichard et al.l (|201Ch . we find that Abell 1835 is 
well described by an unimodal mass distribution. We note 
that the galaxy scale perturbers all present a scale radius 
which is smaller than scale radius inferred for isolated field 
galaxies, in agreement with the tidal stripping scenario. 

The Lenstool software does explore the parameter 
space using a Monte Carlo Markov Chain sampler. At the 
end of the optimization, we have access to these MCMC re- 
alizations from which we can draw statistics and estimate 
error bars. For each realization, we build a two dimensional 
mass map. All these mass maps are then used to compute the 
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Figure 1. Core of Abell 1835. Size of the field equals 82x82 arcsec 2 corresponding to 286x286 kpc 2 . Multiply imaged systems used in 
this work are labeled. 



mean mass map and the corresponding covariance matrix. 
Note that we checked that the PDF of the mass in each pixel 
is approximately Gaussian. These information will then be 
used in the joint fit. 



3 X-RAY DATASETS AND ANALYSIS 

The cluster Abell 1835 is a luminous cluster at redshift 
z = 0.253, which exhibits several indications of a well re- 
laxed dynamical state: for instance its X-ray emission peak 
is associated with a cool core and it is well centred on the 
Brightest Cluster Galaxy (BCG). The X-ray isophotes ap- 
pear quite regular, with a low degree of ellipticity, and with 
the absence of evident substructures. The global (cooling- 
core corrected) temperature T cw has been estimated to be 
T ow = 9.38 ±0.11 keV and an abundance of 0.48 ± 0.03 
solar value ( H3.ip . We classi fy this cluster as a strong cool- 
ing core source fSCC. lMorandi fc Ettorill2007r ). i.e. the cen- 
tral cooling time i CO oi is less than the age of the universe 
tagc.z at the cluster redshift (t C ooi/ia g c,z < 0.1): we estimated 
t coo i ~ 1 x 10 9 yr. As with other SCC sources, Abell 1835 
shows a low central temperature (~ 5 keV) and a strong 



spike of luminosity in the brightness profile. The tempera- 
ture profile is very regular, as expected for relaxed clusters 
(see upper panel of Fig. [2}. 

A full description of the X-ray analysis can be found in 
iMorandi fc Limousin! l|011rj ). Here we only briefly summa- 
rize the most relevant or novel aspects of our data reduction 
and analysis of Abell 1835. 



3.1 X-ray data reduction 

We reduced the Chandra X-ray data using the CI AO data 
analysis package - version 4.3 - and the calibration database 
CALDB 4.4.3. We summarize here briefly the reduction pro- 
cedure. We performed our X-ray analysis on four datasets 
retrieved from the NASA HEASARC archive (observation 
ID 6880, 6881, 7370 and 496) with a total exposure time 
of approximately 200 ks. Three observations (ID 6880, 6881 
and 7370) have been carried out using the ACIS-I CCD 
imaging spectrometer and telemetered in Very Faint mode, 
one (ID 496) using the ACIS-S CCD imaging spectrome- 
ter and telemetered in Faint mode (ID 2321). The level-1 
event files were reprocessed to apply the appropriate gain 
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maps and calibration products and to reduce the ACIS qui- 
escent background. We used the acis_process_events tool 
to check for cosmic-ray background events and to correct 
for eventual spatial gain variations caused by charge trans- 
fer inefficiency to re-compute the event grades. We then fil- 
tered the data to include the standard events grades 0, 2, 
3, 4 and 6 only, and therefore we have filtered for the Good 
Time Intervals (GTIs) supplied, which are contained in the 
f It 1. fits file. The bright point sources were identified and 
masked out using the script vtpdetect, and the result was 
then checked through visual inspection. 

We then used the tool dmextract to create the light 
curve of the background. Indeed a careful screening of the 
background light curve is necessary for a correct background 
subtraction and to discard contaminating flare events. In 
order to clean the datasets of periods of anomalous back- 
ground rates, we used the def lare script, so as to filter out 
the times where the background count rate exceed ±3<r of 
the mean value. Finally, we filtered the ACIS event files on 
energy selecting the range 0.3-12 keV and on CCDs, so as 
to obtain a level-2 event file. 

3.2 X-ray spatial and spectral analysis 

We measure the gas density profile from the surface bright- 
ness recovered by a spatial analysis, and we infer the pro- 
jected temperature profile by analyzing the spectral data. 

The X-ray images were extracted from the level-2 event 
files in the energy range (0.5 — 5.0 keV), corrected by the 
exposure map to remove the vignetting effects, point sources 
were then masked and the images were rebinned by a factor 
of 4 (1 pixel= 1.968 arcsec). 

We determined the centroid (x c ,y c ) of the surface 
brightness by locating the position where the X and Y 
derivatives go to zero, which is usually a robust and outlier- 
resistance approach. We checked that the X-ray emission is 
centered on the BCG: the distance between the X-ray cen- 
troid and the BCG centre is ~ 1.8 arcsec (the uncertainty on 
this measure is comparable to the smoothing scale applied 
to the X-ray image to determine the centroid). 

The spectral analysis was performed by extracting the 
source spectra in circular annuli of radius r* m around the 
X-ray surface brightness centroid. We have selected n* = 8 
annuli out to a maximum distance -R sp cc = 1095 kpc, ac- 
cording to the following criteria: the number of net counts 
of photons from the source in the band used for the spectral 
analysis is at least 2000 per annulus and corresponds to a 
fraction of the total counts always larger than 30 per cent. 
We used the CIAO specextract tool to extract the source 
and background spectra and to construct the redistribution 
matrix files (RMF) and the ancillary response files (ARF). 

Each of the n* annuli hav e been analyzed by using the 
XSPEC package (|Arnaudlll996l ) by simultaneously fitting an 
absorb ed optically-thin plasma emission model (the mekal 
model; [Kaastra 19921; [LTedahl et "ai1ll995l ) to the four obser- 
vations. The fit is performed in the energy range 0.6-7 keV 
by fixing the redshift at z = 0.253, and the photoelectric 
absorption at the galactic value, i.e. to the value inferred 
from radio HI maps. For each of the n* annuli we grouped 
the photons into bins of 20 counts per energy channel and 
applying x 2 -statistics. Thus, for each of the annuli, the free 
parameters in the spectral analysis were the normalization of 



the thermal spectrum K\ oc f dV, the emission-weighted 
temperature Tp ro j i, and the metallicity Zi. 

The four observations were first analyzed individually, 
to assess the consistency of the datasets and to exclude any 
systematic effects that could influence the combined anal- 
ysis. We then proceeded with the joint spectral analysis of 
the four datasets. 

The background spectra have been extracted from re- 
gions of the same exposure for the ACIS-I observations, for 
which we always have some areas free from source emission. 
We also checked for systematic errors due to possible source 
contamination of the background regions. Conversely, for the 
ACIS-S observation we have considered the ACIS-S3 chip 
only and we used the ACIS "blank-sky" background files. 
We have extracted the blank sky spectra from the blank- 
field background data sets provided by the ACIS calibration 
team in the same chip regions as the observed cluster spec- 
tra. The blank-sky observations underwent a reduction pro- 
cedure comparable to the one applied to the cluster data, 
after being reprojected onto the sky according to the obser- 
vation aspect information by using the reproject_events 
tool. We then scaled the blank sky spectrum level to the 
corresponding observational spectrum in the 9-12 keV inter- 
val, since in this band the Chandra effective area is negli- 
gible and thus very little cluster emission is expected. One 
of the advantages of this method is that the derived ARF 
and RMF will be consistent both for the source and the 
background spectrum. However, the background in the X- 
ray soft band can vary both in time and in space, so it is 
important to check whether the background derived by the 
blank-sky datasets is consistent with the real one. From this 
perspective, we verified that for the ACIS-I observations the 
two methods of background subtraction provide very similar 
results for the fit parameters (e.g. the temperature). 



4 SZ DATASET AND ANALYSIS 

The SZ data were collected u sing Bolocam in 200 6, and have 
been presented previously in lSavers et al.l (|201l| y Since that 
publication, these data have been re-reduced using a slightly 
modified reduction pipeline, which we briefly describe here. 
First, the flux calibration model has be en updated based o n 
recent WMAP results as described in ISavers et al.l l|2012f ) , 
which results in < 5% changes to the flux calibration. Sec- 
ond, there were some minor changes to the data-flagging 
procedures, which in general have a very small effect on the 
final SZ image. Third, we note t hat the coord i nates in the 
Abell 1835 image thumbnails in ISavers et all (|201lh were 
offset in declination by 1'9" due to a typo in the source co- 
ordinate file for that cluster. There is no suc h coordinate 
offset for any of the other clusters presented in ISavers et al.l 
(120111) . 

Althou gh the full da t a red uction procedure is described 
in detail in I Savers et al.l (120111 ) , we briefly discuss the rele- 
vant aspects of this processing below. In particular, the data 
are effectively high-pass filtered in a complicated and slightly 
non-linear way in order to subtract noise due to fluctuations 
in the opacity of the atmosphere (i.e., the transfer function 
of the filtering depends weakly on the c luster profile). We 
fit an elliptical gen eralized NFW profile l|Nagai et al. 1 120071 ; 
lArnaud et alj|201oh to Abell 1835, which provided a good fit 
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to the data (x 2 / DOF = 964/945). Note that t his is slightly 
differe nt from the values given in Table 2 of ISavers et al.l 
lj201lf) (x 2 /DOF = 966/945) owing to the slightly differ- 
ent data flagging used in this analysis. The transfer func- 
tion computed for this model was then used for all of our 
subsequent analyses (i.e., all models were filtered using this 
transfer function prior to comparing them to the SZ data). 
We have verified that the biases associated with using this 
single transfer function are negligible compared to the noise 
in the image. The effects of this transfer function can be 
clearly seen in the third panel of Figure [2l where a radial 
profile of the SZ data is plotted. 



5 THREE-DIMENSIONAL STRUCTURE OF 
GALAXY CLUSTERS 

The lensing and the X-ray/SZ emission both depend on the 
properties of the DM gravitational potential well, the for- 
mer being a direct probe of the two-dimensional mass map 
via the lensing equation and the latter an indirect proxy of 
the three-dimensional mass profile through the HE equation 
applied to the gas temperature and density. In order to infer 
the model parameters of both the IC gas and of the underly- 
ing DM density profile, we perform a joint analysis of SL and 
X-ray/SZ data. We briefly outline the methodology in order 
to infer physical properties in triaxial galaxy clusters: (f ) We 
start with a generalized Navarro, Frenk and White (gNFW ) 
triaxial model of the DM as described in iJing fe Sutd(|2002h . 
which is representative of the total underlying mass distri- 
bution and depends on a few parameters to be determined, 
namely the concentration parameter C, the scale radius P s , 
the inner slope of the DM 7 , the two axis ratios (rjnM,a 
and r/pM.b) and t h e Eul er angles ip, 6 and <f> (2) following 
iLee fc Sutol l|2003l . I2004T ). we recover the gravitational po- 
tential and two-dimensional surface mass S (Equation If f I) 
of a dark matter halo using this triaxial density profile; (3) 
we solve the generalized HE equation, i.e. including the non- 
thermal pressure P nt (Equation [SJ, for the density of the IC 
gas sitting in the gravitational potential well previously cal- 
culated, in order to infer the theoretical three-dimensional 
temperature profile T; (4) we calculate the SZ temperature 
decrement map AT(v) (Equation [9} and the surface bright- 
ness map Sx (Equation |SJ) related to the triaxial ICM halo; 
and (5) the joint comparison of T with the observed temper- 
ature, of Sx with the observed brightness image, of AT(u) 
with the observed SZ temperature decrement, and of £ with 
the observed two-dimensional mass map gives us the param- 
eters of the triaxial ICM and DM density model. 

Here we briefly summ arize the major findings of 
iMorandi fc Limousin! (|01ld ) for the triaxial joint analysis 
in order to infer the desired physical properties, as well as 
the improvements add ed in the curren t analysis; additio nal 
details can be found in IMorandi et all (|2007l , |2010| , lOllal lbh ; 
IMorandi fc Limousin] (lOl id). 

In Morandi fc Limousin! (|01ld ) we model the DM and 
ICM ellipsoids to be orientated in an arbitrary direction on 
the sky. We introduced two Cartesian coordinate systems, 
x = (x,y,z) and x' = (x',y',z'), which represent respec- 
tively the principal coordinate system of the triaxial dark 
halo and the observer's coordinate system, with the origins 
set at the center of the halo. We assumed that the «'-axis lies 



along the line of sight to the observer and that the x' , y' axes 
identify the direction of West and North, respectively, on the 
plane of the sky. We also assumed that the x, y, z-axes lie 
along the minor, intermediate and major axis, respectively, 
of the DM halo. If we define tj), 6 and cj> as the rotation 
angles about the x, y and z axis, respectively, then the rela- 
tion between the two coordinate systems can be expressed 
in terms of the rotation matrices M x (ip), M y (9), M z ((j>) with 
Euler angles ip,6,cj): 



c' = M x {il>)#M y {e)#M z {ct>)#y 



(1) 



In order to parameterize the cluster mass distribution, 
we consider a triaxial generalized Navarro, Frenk & White 
model gNFW (IJing fc Sutdl2002D : 



p(R) = 



ScPC,: 



(R/Rs) 1 (1 + R/Rs) 



3 — 7 : 



(2) 



where R s is the scale radius, 5c is the dimensionless char- 
acteristic density contrast with respect to the critical den- 
sity of the Universe pc,z at the redshift z of the cluster, 
and 7 represents the inner slope of the density profile; 
pc,z = 3H(z) 2 /8tvG is the critical density of the universe 



at redshift z, H z = E z H Q , E z = [Q M (1 + z) 3 + A 
200 C 3 



1/2 



&c 



3 F(C, 7 ) 



and 



(3) 



where C = R200/RS is th e concentration para meter, and 
F(C,y) has been defined in lWvithe et~aH l|200lh . 

The radius R can be regarded as the major axis length 
of the iso-density surfaces: 



y_ 

b 2 



(a ^ b ^ c). 



(4) 



We also define r\oM,a = a/c and T]DM,b = b/c as the minor- 
major and intermediate-major axis ratios of the DM halo, 
respectively, and et and e c the relative eccentricities (e.g. 

e b = y/1 -(b/c) 2 ) 

The work of ILee fc Sutol (|2003r ) showed that the iso- 
potential surfaces of the triaxial dark halo are well approx- 
imated by a sequence of concentric triaxial distributions of 
radius Pi C m with different eccentricity ratio. For Pi cm it 
holds a similar definition as R (Equation 0]) , but with IC 
gas eccentricities e b and e c . Note that e b = eb(eb,u, 7) and 
e c = £c(e , u, 7), with u — R/R s , unlike the constant eb,e c 
for the adopted DM halo profile. In the whole range of it, 
f-b/e-b (fic/ec) is less than unity (~ 0.7 at the center), i.e., 
the intracluster gas is al together more s pheric al than the 
underlying DM halo (see IMorandi et al.l (|20f 01 ) for further 
details) . 



5.1 X-ray, SZ and lensing equations 

For the X-ra y analysis we rely on a generalization of the 
HE equation (|Morandi et al.ll"011bh . which accounts for the 
non-thermal pressure P nt and reads: 



VPtot = -PgasV$ 



(5) 



where p gas is the gas mass density, $ is the gravitational 
potential, P tot = P t h + Pnt ■ We implemented a model where 
Pnt is a fraction of the total pressure P to t, and we set this 
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fraction to be a power law with radius (|Shaw et alj feoiO): 



= £ (R/Saoo) 

-Ttot 



(6) 



Note that X-ray and SZ data probe only the thermal com- 
ponent of the gas P t h = n e kT, k being the Boltzmann 
constant. From Equations ([5]) and (|6]) we point out that 
neglecting P nt (i.e. P tot = P t h) systematically biases low 
the determination of cluster mass profiles. We stress that 
th e model in Equation (l6[) is an improvement with respect 
to lMorandi &: Limousinl i Olid ), where we assumed that the 
non-thermal pressure is a constant fraction of the total pres- 
sure. 

Given that Equation ([5]) is a first order differential equa- 
tion, we need a boundary condition on the pressure, P, 
which represents the pressure at -R200, and it is an unknown 
parameter to be determined. 

To model the electron density profile in the triaxial ICM 
halo, we use the following fitting function: 

n e (R icm ) = n (R icln /r ci y s (l + Rf cm /r 2 cl )- 3/2e+S/2 ■ 

■(l + fltan/rt)~" /4 (7) 

with parameters (no, r ci , e, 8, r C2 , v). Note that the fitting 
function in J7]) has more degree s of fre edom than that em- 
ployed in iMorandi fc Limousinl |011cl ). We computed the 
theoretical three-dimensional temperature T by numerically 
integrating the equation of the HE (Equation [5]), assuming 
triaxial geometry and a functional form of the gas density 
given by Equation (J7J). 

The observed X-Ray surface brightness Sx is given by: 



Sx = 



47r(l + z) 



-A(T p * roj , Z) I n c n p dz' , 



(8) 



where A(T* I0 j, Z) is the cooling function. Since the projec- 
tion on the sky of the plasma emissivity gives the X-ray sur- 
face brightness, the latter can be geometrically fitted with 
the model n e (-Ricm) of the assumed distribution of the elec- 
tron density (Equation [7|) by applying Equation {SJ. This 
has been accomplished via fake Chandra spectra, where the 
current model is folded through response curves (ARF and 
RMF) and then added to a background file, and with absorp- 
tion, temperature and metallicity measured in the neighbor- 
ing ring in the spectral analysis ( H3.2p . In or der to calculate 
A(r* rni ,Z), we ad opted a MEKAL model (jKaastral [l99l 
iLiedahl et alll 19951 ) for the emissivity. 

The thermal SZ effect is expressed as a small variation 
in the temperature AT(v) of the CMB as a function of the 
observation frequency: 



ATM 



J P{v)f{u-T{v))dz' 



(9) 



where ot is the Thomson cross-section, P e (r) = 
n e (r)fc{,T e (r) is the pressure of the electrons of the ICM at 
the volume element of coordinate r, kb is the Boltzmann 
constant, I = 2(kT cmb )' i / (he) 2 , T cmb = 2.725 K. 

f(v;T(r)) takes into account the spectral shape of the 
SZ effect and it reads: 



/(^;T(r)) = (x 



+ 1 



•4)*(l + 0/ (x;T)) 



(10) 



where x = hu/kT cm ^ accounts for the frequency dependence 
of the SZ effect, and for the relativistic corrections related to 



the term Of(x,T) |ltoh et al.lll998l ). Note that in Equation 
([9]) we account for the implicit dependence of f(v;T(r)) on 
radius. 

Next, the two-dimensional SZ model AT(u) is con- 
volved with the Bolocam point-spread function and the mea- 
sured transfer function. In practice, the transfer function 
convolution is performed via multiplication in the Fourier 
domain. This filtering significantly reduces the peak decre- 
ment of the cluster and creates a ring of positive flux at 
r ~ 2 arcmin. This filtered model is then compared to the 
observed SZ temperature decrement map. We also calcu- 
lated the noise covariance matrix C among all the pixels of 
the observed SZ temperature decrement map through 1000 
jackknife realizations of our cluster noise. In this perspective 
we assumed that the noise covariance matrix for the SZ data 
is diagonal, a s this was shown to be a good assumption in 
ISavers et all (|201ll ). 

For the lensing analysis the two-dimensional surface 
mass density £ can be expressed as: 



p{R)dz' 



(11) 



We also calculated the covariance matrix C among all the 
pixels of the reconstructed surface mass fsee lMorandi et aU 
|011bl ) for further details). 



5.2 Joint X-ray+SZ+lensing analysis 

The probability distribution function of model parame- 
ters has been evaluated via a Markov Chain Monte Carlo 
(MCMC) algorithm, by using the likelihood £ described be- 
low as proposal density and a standard method for rejecting 
proposed moves. This allows to compare observations and 
predictions, and to infer the desired physical parameters. 
The likelihood has been constructed by performing a joint 
analysis for SL and X-ray/SZ data. More specifically, the 
system of equations we simultaneously rely on in our joint 
X-ray+SZ+Lensing analysis is: 



T(C, Rs, 7, l]DM,a, TjDM,b, 1p, 0, (f>, 7l , f cl , £, S, T C2 , V, £, 71, P) 

Sx(C, R B ,j, T]DM,at TjDM,bi ip,0,<t>, no,r cl , e, 5, r C2 , v) 
AT(C, -R s , 7, r}DM,a,VDM,b, ip, 6, 4>, n , r ci , s, S, r C2 , v, £, n, P) 
S(C, Rs, 7, r) DM ,a, rjDM.b, ip, 8, <f) (12) 

where the parameters C (concentration parameter), R s 
(scale radius), 7 (inner DM slope), r/DM,a (minor- major axis 
ratio), r]DM,b (intermediate-major axis ratio), and ip,9,(j> 
(Euler angles) refer to the triaxial DM halo (Equation [5}; 
the parameters no, r ci , e, S, r C2 , v refer to the IC gas density 
(Equation [7|); £,n (normalization and slope, respectively) 
refer to the non-thermal pressure (Equation O; and P to 
the pressure at -R200, which is a boundary condition of the 
generalized HE equation (Equation [5]). 

In our triaxial joint analysis the three-dimensional 
model temperature T is recovered by solving equation 
and constrained by the observed temperature profile, the 
surface brightness is recovered via projection of the gas den- 
sity model (Equation [S} and constrained by the observed 
brightness, the SZ signal is deduced via projection of the 
three-dimensional pressure (Equation [9]) and constrained by 
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the observed SZ temperature decrement, and the model two- 
dimensional mass S is recovered via Equation and con- 
strained by the observed surface mass. 

Hence the likelihood L oc exp(— x 2 /2), and x 2 reads: 

X = Xx,T + Xx,S + Xsz + Xkms (13) 

with Xx,t. Xx,s. Xsz and xL s equal to th e X 2 coming from 
the X-ray temperature, X-ray brightness, SZ temperature 
decrement and lensing data, respectively. 

For the spectral analysis, x\ t is equal to: 

Xx.t = £ V (14) 

1—1 proja 

^proj.i being the observed projected temperature profile in 
the ith circular rin g and Tproi.i the azim uthally-averaged 
projection (following M azzotta et al.ll2004h of the theoretical 
three-dimensional temperature T; the latter is the result of 
solving the HE equation, with the gas density n e (J?i cm ). 
For the X-ray brightness, Xx,s reads: 

2 sr^x^ (Sx,i — Sx,i) 2 , , 

xx,s = LL -2 ( 15 ) 

j 1=1 a s,i 

with Sx,i and S x ,i equal to the theoretical and observed 
counts in the ith pixel of the jth image. Given that the 
number of counts in each bin might be small (< 5), then 
we cannot assume that the Poisson distribution from which 
the counts are sampled has a nearly Gaussian shape. The 
standard deviation (i.e., the square-root of the variance ) for 
this low-count case has been derived bv lGehrelsl (|l986h : 

o s ,i = 1 + y/S' Xii + 0.75 (16) 

which has been demonstrated to be accurate to approxi- 
mately one percent. Note that we added background to Sx,i 
as measured locally in the brightness images, and that vi- 
gnetting has been removed in the observed brightness im- 
ages. 

For the SZ (lensing) constraint D, the Xd contribution 

is: 

X 2 D = [D-D*] t C 1 [D-D*] , (17) 

where C is the covariance matrix of the two-dimensional SZ 
temperature decrement (projected mass), D* are the ob- 
served measurements of the two-dimensional SZ tempera- 
ture decrement (projected mass) in the ith pixel, and D is 
the theoretical 2D model. Note that we removed the central 
25 kpc of the 2D projected mass in the joint analysis, to 
avoid the contamination from the cD galaxy mass. 

We report the average value and standard deviation of 
the marginal probability distribution for the individual pa- 
rameters. In addition to a complete statistical analysis of 
the chain, we performed a ser ies of convergence tests : the 
Gelman & Rubin R statistics (|Gelman fc Rubm]|l992l ) and 
a splittest (which essentially consists in splitting the chain 
into 2, 3 or 4 parts and comparing the difference in the 
parameter quantiles). We confirmed the convergence of our 
result using these tests. 

We can determine the physical parameters of the clus- 
ter, for example the 3D temperature T, or the shape of the 
DM and the ICM, by relying on the generalized HE equa- 
tion and on the robust results of the hydrodynamical sim- 
ulations of the DM profiles (i.e. gNFW). We also point out 



that, given the complementary data sets that have been in- 
cluded in this analysis, we do not need to rely on any prior 
from theoretical predictions, such as priors on the concen- 
tration parameter, on the halo mass or on the axis ratio (e.g. 
ICorless et alj|2009f ). which might be biased due our incom- 
plete understanding of the cluster physics in simulations. 

In Fig. [2] we present an example of a joint analysis for T, 
Sx, AT(y) and E: for Sx, AT(u) and E the ID profile has 
been presented only for visualization purposes, the fit being 
applied on the 2D X-ray brightness/SZ/surface mass data. 
Note that in the joint analysis the X-ray, SZ and lensing 
data are all well described by our model, with a Xrcd = 
1.04(1477928 degrees of freedom). 



6 RESULTS AND DISCUSSION 

In the previous section we showed how we can determine the 
physical parameters of the cluster by fitting the available 
data relying on the HE equation and on a DM model that 
is based on robust results of hydrodynamical cluster simula- 
tions. In this section we present our results and discuss their 
main implications. We particularly focus on the implications 
of our analysis for the determination of the full triaxiality, 
viability of the CDM scenario, the presence of non-thermal 
pressure and the gas properties in the outskirts. 

6.1 Model parameters 

The model parameters are summarized in table [5] Our 
work indicates that Abell 1835 is a triaxial galaxy clus- 
ter with DM halo axial ratios r]oM,a ~ 0.59 ± 0.05 and 
VDM,b = 0.71 ± 0.08, and with the major axis slightly in- 
clined with respect to the line of sight of 8 — 18.3 ± 5.2 
deg. Our findings strengthen the view of a triaxial cluster 
elongated along the line of sigh t , in a greement with the pre- 
dictions of lOguri fc Blandfordl l|2009h . who showed that SL 
clusters with the largest Einstein radii constitute a highly bi- 
ased population with major axes preferentially aligned with 
the line of sight, thus increasing the magnitude of the lensing 
signal. 

The axial ratio of the gas is r? g as,a ~ 0.77—0.86 and 
7?gas,b ~ 0.79—0.87, moving from the center toward the virial 
radius. 

The value of the virial radius is: -R200 = 3809 ± 254 
kpc. Note that we used a 'triaxial' definition of R200, which 
refers to the major axis (roughly along the line of sight) of 
the triaxial DM halo (see Equation^, so to make a compar- 
ison with a scale length on the plane of the sky we should 
multiply -R200 by ~ r]DM,a, i.e. the virial radius corresponds 
to a scale length of ~ T]DM,a ■ R200 ~ 2240 kpc. This also 
means that, given the SZ measurements out to ~ 2200 Kpc, 
SZ constrains the IC gas out to ~ -ffeoo- 

Another main result of our work is the estimate of the 
non-thermal pressure support, at a level up to ~20% in the 
outer volumes (~ ifeoo)- 

In order to assess the importance of SZ data, we per- 
formed the following test. We excluded the SZ data from 
the joint analysis, fixing P to the pressure at the X-ray 
boundary and by assumin g the model of n e employed in 
iMorandi fc Limousin! (|011cT ): this different modelling of n e is 
needed because the parameters r C2 , v and n are very poorly 
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Figure 2. Example of the joint analysis of T, Sx, AT(i/) and X. In the upper panel we display the two quantities which enter in 
the X-ray analysis (Equation 1 14jl : the observed spectral projected temperature Tp ro j m (big points with errorbars) and the theoretical 
projected temperature T pro j m (diamonds). We also show the theoretical 3D temperature (solid line), which generates T pro j m through 
convenient projection techniques. In the second panel from the top we display the two quantities which enter in the X-ray brightness 
analysis (Equation ll5ll : the observed surface brightness profile S x (points with errorbars) and the theoretical one Sx (solid line). In the 
third panel from the top we display the two quantities which enter in the SZ temperature decrement analysis (Equation [9{: the observed 
SZ temperature decrement profile (points with errorbars) and the theoretical one AT(i/) (solid line). Both the observed and theoretical 
SZ temperature decrement are convolved with the transfer function: note that this filtering significantly reduces the peak decrement of 
the cluster and creates a ring of positive flux at r ~ 2 arcmin. In the lowest panel we display the two quantities which enter in the 
lensing analysis (Equation 1171 ): the observed surface mass profile S* (points with error bars) and the theoretical one X (solid line). 
Note that for surface brightness (surface mass) and the SZ data the ID profile has been presented only for visualization purposes, the 
fit being applied on the 2D data. Moreover, for the surface brightness we plotted data referring to the observation ID 6880. The virial 
radius corresponds to a scale length on the plane of the sky of ~ r]r)M a ' -R200 ~ 2240 kpc. 
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Table 2. Model parameters of Abcll 1835. Lines 1 — 8 refer to 
the model parameters C (concentration parameter), R B (scale ra- 
dius), 7 (inner DM slope), "cm, a (minor-major axis ratio), r)pM,b 
(intermediate-major axis ratio), and ip,8,<j> (Euler angles) of the 
DM halo (Equation [2). The lines 9 — 14 refer to the model param- 
eters no, r c . , e, 8, r C2 , v of the IC gas density (Equation 0, while 
the lines 15 — 16 to the model parameters £, n (normalization 
and slope, respectively) of the non-thermal pressure (Equation 
Finally, the last line refers to the model parameter P of the 
pressure at i?200 1 which is a boundary condition of the generalized 
HE equation (Equation . 



c 


4.32 ± 0.44 


Rs (kpc) 


891.0 ± 114.3 


7 


1.01 ±0.06 


VDM.a 


0.59 ±0.05 


VDM,b 


0.71 ±0.08 


i> (deg) 


-55.0 ±6.9 


6 (deg) 


18.3 ±5.2 


4> (deg) 


3.8 ±4.6 


n (cm -3 ) 


0.018 ±0.002 


r ci (kpc) 


117.7 ±10.1 


e 


0.68 ±0.02 


8 


0.82 ±0.03 


r C2 (kpc) 


1674.3 ± 266.7 


V 


0.44 ±0.04 


e 


0.177 ±0.065 


n 


0.77 ±0.21 


P (erg/cm 3 ) 


(2.7 ±0.7) x 10- 



constrained without the SZ datasets. We obtained larger 
(10%-25%) the errors on the final parameters with respect 
to the case where we include SZ data. From this test, it is 
clear how much SZ data are important to remove degen- 
eracy among the parameters and crucial in measuring the 
physical properties of the IC gas in the outskirts, which are 
inaccessible to X-ray and SL observations. Therefore, in the 
present analysis we need to jointly combine all the datasets 
(X-ray, SZ and lensing) to determine the desired physical 
parameters. 

We stress that in the internal regions the physical prop- 
erties of the cluster are overconstrained by our data (e.g., the 
thermal pressure is measured directly from the SZ data and 
also by the combination of X-ray density and spectroscopic 
temperature). This provides critical insights to our under- 
standing of clusters, and critical tests of current models for 
the formation and evolution of galaxy clusters. Yet, in the 
outskirts we mainly rely on the SZ data, since the outer vol- 
umes are not constrained by X-ray and SL observations. As 
a note of caution, we point out that the physical properties 
in the outskirts are then no longer overconstrained. While 
SZ data are well described by our model, adding other con- 
straints, e.g. weak lensing (WL) data and/or Suzaku X-ray 
data, might help in gauging the impact of potential system- 
atics on the desired properties in the outer volumes. 

In Figure [3] we present the joint probability distribu- 
tion among different parameters in our triaxial model. For 
example, we point out that there is a positive correlation 
between £ and ??DM,a, i.e. the X-ray/Lensing mass discrep- 
ancy in clusters with prominent strong lensing features can 



be explained via a combination of both triaxiality and non- 
thermal support. 

We also tried to gauge the impact of possible system- 
atics on the inferred physical parameters. For example, the 
choice of using the dPIE mass distribution for lensing-only 
data (which constraints the 2D mass out to ~ 300 kpc) 
might affect, in principle, the derived parameters, for exam- 
ple the non-thermal component (which is constrained by the 
data out to ~ J?20o). In order to test this assumption, we 
also fit the lensing-only data with a gNFW model, and we 
found a projected mass profile consistent with the one de- 
rived from a dPIE mass distribution within a few percent. 
Actually, the small range constrained by the SL data does 
not allow us to discriminate between a dPIE and a gNFW 
profile. This means that the actual systematic uncertainty 
on the physical parameters does not depend on the assumed 
model of DM for lensing-only data (dPIE or gNFW), but 
mostly stems from calibrations uncertainties of X-ray and 
SZ data. For X-ray the calibration un certainties have been 
estimated at ~6% l|Reese et al.ll2010j), while for Bolocam 
Sunyaev Zeldovich data in ~5% (|Savers et alJlioTH . l2012h . 
We then repeated the MCMC analysis by including the pre- 
vious systematics in the X-ray and SZ data (assuming that 
they have Gaussian distributions): we found that the errors 
on the inferred parameters get slightly larger (~ 10 — 20 
percent). 

Finally we also compared the azimuthal angle <f> = 
76.4 ± 4.7 deg and the eccentricity on the plane of the 
sky (e = 0.16 ± 0.03), with the values on the total two- 
dimensional mass from the analysis from Lenstool (4> = 
77.4 ± 0.6 deg and e = 0.11 ± 0.02): note the good agree- 
ment. Fo r the method to reco ver e and <j> we remand the 
reader to iMorandi et "ail i|2010T ). 

6.2 Implications for the viability of the CDM 
scenario 

Clusters are an optimal place to test the predictions of cos- 
mological simulations regarding the mass profile of dark 
matter halos. A central prediction arising from simulations 
of cosmic structure formation in a hierarchical, dark matter 
dominated universe is that the density profile of DM halos 
is universal across a wide rang e of mass scales fro m dwarf 
galaxies to clusters of galaxies l|Navarro et al.lll997| ): within 
a scale radius, r s , the DM density asymptotes to a shallow 
powerlaw trend, pdm{t) oc r~ 7 , with 7 = 1, steepening at 
increasing radii. Simulations also suggest that galaxy clus- 
ter concentrations, which are a measure of a halo's central 
density, decrease gradually with virial mass. Nevertheless, 
the value of the logarithmic inner slope 7 and the actual 
mass-concentration relation are still debated. 

Recent works investigating mass distributions of in- 
dividual galaxy clusters have measured high concentra- 
tion parameters, very large Einstein radii, and high effi- 
ciency in generating giant arcs, which represents a ma- 
j or inconsistency with t he theoretical CDM expectations 
dBroadhurst et ail [2005: Lim ousin et al. I l007al : IZitrin et al.l 
iQllal). 

Moreover, measurements of 7 over various mass scales 
have proved controversial, yielding conflicting values of 7, 
with large scatter from one cluster t o another. Shallow cusps 
have been inferred by the analysis of lSand et al.l |2008) , rais- 
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Figure 3. Marginal probability distribution among different parameters in our triaxial model. The solid(dashcd) line represent the l(2)-a 
error region, while the big point represents the expectation value. 
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ing doubts on the predictions of the CDM scenario, but these 
determination rely on the standard spherical modeling of 
galaxy clusters. 

The disagreement between theory and observation 
might be explained by triaxiality. In clusters that are elon- 
gated along the line of sight the measured concentration 
parameter is substantially biased up with respect to the 
theoretical expect ations, and the observe d lensing proper- 
ties are boosted |Meneghetti et al-llioloh . Moreover, elon- 
gation/flattening of the sources along the line of sight, as 
well as the degeneracy of 7 with other parameters, i.e. the 
concentration parameter and the scale radius, likely affect 
the estimated values of 7. Therefore knowledge of the in- 
trinsic shape and orientation of halos is crucial to unbiased 
determinations of inner slope of the DM and concentration 
parameter, and hence to assess the viability o f the standard 
cosmological framework (|Morandi et al-lfeoiOl ). 

One of the main results of the work presented here is a 
measurement of a central slope of the DM 7 = 1.01±0.06 by 
means of a triaxial joint analysis for Abell 183 5: this value is 
in agr eement with the CDM predictions from lNavarro et"ahl 
l|l997h (i.e. 7 = 1). We point out that we removed the cen- 
tral 25 kpc in the lensing data, to avoid the contamination 
from the cD galaxy, though we checked that there is a very 
little dependence of the physical parameters on the radius of 
the masked region. The value of the concentration param- 
eter C — 4.32 ± 0.44 is in agreement with the theoretica l 
expectation from N-body simulations of iNeto et al.l (|2007h : 
iDuffv et all (|2008T ). where C ~ 4 at the redshift and for the 
virial mass of Abell 1835, and with an intrinsic scatter of 
~ 20 per cent. By means of a lensing-only analysis under 
the assumption of spherical geometry, we infer a large value 
of the concentration C = 6.19 ± 0.63, which lies above the 
standard C—M relation, putting some tension with the pre- 
dictions of the standard model. This confirms our insights 
about the role of the effects of geometry on the cluster con- 
centrations. 

Given that numerical simulations customarily retrieve 
the concentration parameter by using a spherically-averaged 
total density profile, the question arises whether a compar- 
ison of such simulations with the concentration parameter 
inferred in our triaxial framework is tenable. In this per- 
spective, we generated an ideal ellipsoidal NFW cluster with 
parameters fixed to those of Abell 1835 (we fixed 7 = 1 for 
simplicity in the comparison with simulations), and then we 
measured C by spherically-averaging the total density pro- 
file. We found that such concentration parameter is slightly 
lower (~ 4%) with respect to the value in the triaxial frame- 
work. Given that this bias is much smaller than t he intrinsic 
scatt er of C in numerical simulations (~ 20%, see lNeto et al.l 
120071 ). we conclude that a comparison of our findings with 
simulations is still convincing. 

We also report the value of the concentration parameter 
for an X-ray-only analysis under the assumption of spherical 
geo metry: C = 4.40 ± 0.23 . 

IMorandi et alj (|2010l . lOllal lbh; iMorandi fc Limousin! 
(|01 lq ) analyzed the galaxy clusters Abell 1689 and 
MACS J1423. 8+2404 in our triaxial framework: we found 
that 7 and C are close to the CDM predictions for these 
clusters, in agreement with the results in the present pa- 
per. Our findings provide further evidences that support the 
CDM scenario. 



We also find a very good agreement of our inferred 
value of the concen tration parameter and axial ratios with 
ICorless et al.l (|2009h . who constrained the triaxial shape of 
the total mass distribution of Abell 1835 via weak lensing 
data and under a range of Bayesian priors derived from the- 
ory. 

6.3 Non-thermal gas pressure 

Hydrodynamical simulations have shown that a significant 
fraction of the total energy of the IC gas is non-thermal, 
mainly due to random gas motions and turbulence in the 
same IC gas. This energy is sourced by several mecha- 
nisms: plasma instabilities, mergers and subcluster assem- 
bly in the hierarchical structure formation scenario, shock 
waves, wakes of galaxies moving into the IGM, outflows 
from active galactic nuclei (AGNs ) hosted in the center o f 
galaxy clusters and galactic winds 1 Norman Sz Brvanlll999l '). 
In particular, in the hierarchical structure formation sce- 
nario turbulent motions should occur in the IC gas while 
the matter continues to accrete along filaments. Indeed gas 
accreting onto clusters of galaxies has bulk velocities of 
abou t v = 1900(T/6.7) ' 52 km/s at 1 Mpc (see lMiniati et all 
l2000t ). T being the gas temperature. This velocity is com- 
parable to the expected sound speed of 1000-1500 km/s of 
the ICM, and hence generates turbulent gas motions at the 
boundary between the bulk flow and the thermalized ICM 
l|Vazza et al.ll2009l ) . This non- hydrostatic energy should then 
cascade from large to small scales and can eventually dissi- 
pate into the gas, leading to a (partial) thermalization of 
the IC gas. Yet, the total energy budget in the form of tur- 
bulent motions inside galaxy clusters, as well as their dis- 
tribution and their connection with cluster dynamics and 
non-gravitational processes in galaxy clusters is still open to 
debate in the literature. 

While the level of non-thermal pressure is typically 
found to be small in the central regions of clusters, it in- 
creases with radius, becoming a sig nificant fraction of the to- 
tal pressure in the outer volumes l|Lau et al.ll2009i ). It is also 
clear that non-thermal pressure support causes a systematic 
underestimate of the c luster mass recovered under the as- 
sump tion of strict HE |Nagai et al.ll2007l : iMeneghetti et al.l 
2010). This translates into biases in the determination of the 
cosmological parameters. 

A number of observational evidences of non-thermal 
pressure have also be en published in the last few years. 
ISchuecker et~afl (|2004 ) obtained spatially-resolved gas pres- 
sure maps of the Coma cluster which indicated the pres- 
ence of a significant amount of turbulence, with a spectrum 
of the fluctuations consistent with Kolmogorov turbulence. 
This yielded the lower limit of ~10% of the total IC gas 
pressure in turbulent form. Additional evidences of turbu- 
lent motions inside nearby galaxies come fro m the obser- 
vatio n of pressure fluctuations in Abe ll 754 (|Henrv et al.l 
l2004h and Perseus dFabian et all 120031 ) . The observational 
results of iMahdavi et al. ( 20081 ) based on X-ray and WL 
mass determinations indicate that there is a radial trend of 
the X-ray /WL mass ratio, that is interpreted as caused by 
non-thermality increasing toward the outer regions. Never- 
theless, their findings hinge on the assumed spherical geom- 
etry, so they did not disentangl e the effect of triaxiality from 
non-thermal pressure support. IMorandi fc Limousinl (|01ld ) 
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measured the non-thermal component of the gas (~ 10%) 
relative to the total energy budget of the inner volumes of 
the IC gas of Abell 383. 

Motivated by the need to study the magnitude and gra- 
dient of the non-thermal pressure support and given the con- 
straints on the IC gas provided by SZ data out to ~ -R200, 
here we implemented a model where P n t is a fraction of the 
total pressure P to t , and we set this fractio n to be a power law 
with the radius. The theoretical work of IShaw et al.l (|2010l ) 
based on 16 simulated clusters shows that this model reason- 
ably reproduces the trend of P n t/ftot throughout a cluster. 
They found a slope of 0.8±0.25, in agreement with our mea- 
sured value, and a normalizatiorjf] > 0.3, in marginal (~ 2a) 
disagreement with our analysis, where we found a normaliza- 
tion £ = 0.177 ± 0.065. Therefore our findings indicate that 
the level of turbulence in the numerical simulations might 
be overestimated. 

We stress that this is the first observational measure- 
ment of the non-thermal pressure out to R200 , recovered via 
a triaxial joint X-ray, SZ and lensing analysis. Our results 
therefore can provide an anchor for numerical models of ICM 
physics and for simulations of the formation and ongoing 
growth of galaxy clusters, given that measurements of the 
non-thermal energy in the IC gas are a proxy of the amount 
of energy injected into clusters from mergers, hierarchical 
assembly of substructures, accretion of material or feed- 
back from AG Ns. In this per s pectiv e, the cooling-preheating 
simulations of IStanek et al.l (|2010f) suggest that the IC gas 
in the gravitation-only simulations develops more substruc- 
tures with time than in the former. The suppression of sub- 
structures caused by the preheating leads to a lower level of 
kinetic energy in bulk motions with respect to simulations 
without preheating, confirming that the level of non-thermal 
pressure is sensitive to the particulars of the physical pro- 
cesses. Therefore a more extensive physical treatments that 
incorporates further physical processes in the ICM which 
are currently uncertain, for example galaxy and supermas- 
sive black hole formation, MHD, viscosity, conduction, star 
formation feedback, magnetic fields and non-thermal plas- 
mas, (pre)-heating, might be needed (or improved) in sim- 
ulations to match the amount of non-thermal energy with 
our observational findings. 

Moreover, given that hydrodynamical simulations indi- 
cate that non-thermal pressure provided by gas motions sig- 
nificantly modifi es the ICM profiles in the cluster outskirts 
l|Lau et alj|2009l ). the non-thermal pressure must be accu- 
rately determined in order to obtain unbiased measurements 
of the physical parameters. 

As a note of caution, we remind the reader that our 
results were obtained for just one galaxy cluster, though it 
is expected that the observed physical processes are common 
at least for relaxed objects. In this perspective we plan to 
collect data for a larger sample of clusters to strengthen the 
statistical significance of our fin dings. 

We also report the work of ISanders et al.l ((2010) , who 
placed a direct limit on turbulence based on the non-thermal 
velocity broadening measured from the emission lines origi- 
nating in the central 30 kpc of Abell 1835. They found that 



3 Note that in lShaw et al.l d2010f) the scale radius for P n t refers 
to i?500i so we renormalize their results to -R200- 
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Figure 4. Top panel: spherically-averaged entropy profiles S for 
Abell 1835 recovered via triaxial joint X-ray, SZ and lensing anal- 
ysis. The solid line represents the expectation value for S, while 
the 1-cr errors are represented by the gray shaded regi on. The 
dashed line represents the predictions of lVoit et~a l. (005b), where 
the entropy is defined as S(r) = S200 l-32(r/r2oo) 11 i 5200 being 
a characteris tic value of the en tropy at an overdensity of 200 (see, 
e.g., Eq. 2 in lVoit et alj|005bl) . Lower panel: spherically-averaged 
gas density profile n e for Abell 1835. The solid line represents 
the expectation value for n e , while the 1-cr errors are represented 
by the gray sha ded region. The long d ashed line represents the 
predictions from Roncarelli et al. (2006). 



the ratio of turbulent to thermal energy density in the core 
is less than 13%, in agreement with the present work. 



6.4 Physical properties in the outskirts 

The outskirts of galaxy clusters present an opportunity to 
study the formation of large scale structure as it happens. 
They have special importance in cluster cosmology, because 
they are believed to be much less prone to complicated clus- 
ter astrophysics, such as radiative gas cooling, star forma- 
tion, and energy injection from AGNs, although they are po- 
tentially more susceptible to the turbulence and bulk flows 
that result from structure formation processes (3H3I- The 
physics of the IC gas in the outer volumes is relatively sim- 
ple and nearly self-similar, being dominated by the gravity- 
driven collisionless dynamics of DM and hydrodynamics of 
the gas. However, until very recently, observational studies 
of the ICM have been limited to radii considerably smaller 
than the virial radius of clusters. Here we aim at understand- 
ing the properties of the ICM in the outskirts of Abell 1835, 
comparing our findings with the results of hydrodynamical 
numerical simulations. 

Our spherically-averaged gas density has good agree- 
ment with the predictions from hydrodynamical numeri- 
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cal simulations inc luding cooling, star fo rmation and su- 
pernovae feedback (jRoncarelli et al.ll2006l ). although with a 
slightly flatter slope. A possible explanation of this trend 
is that in the above simulation Roncarelli and collabora- 
tors did not include AGN feedback and/or preheating, which 
might be important even at large radii, smoothing the ac- 
cretion pattern and leading to a flatter gas distribution 
l|Borgani et al.ll2005l ). 

It is interesting to point out that some recent Suzaku ob- 
servations indicate shallow density /entropy profiles in clus- 
ter outskirts, at varian ce with the re s ults fr om previous 
ROS AT observations l|Vikhlinin et all QUI; lEckert et all 
I2011T), and with the results from numerical simulations 
1 Roncarelli et alj|2006j ). Thus, the behavior of the gas den- 
sity in clust er outskirts is still t he su bj ect of debate. In this 
perspective, ISimionescu et al.l (120111 ); iGeorge et al.l (2009) 
analyzed Suzaku X-ray observations along narrow arms, and 
they found that the electron density decreases steadily with 
radius, approximately following a power-law model, while we 
observe a general trend of steepening in the radial profiles 
of the gas density beyond 0.2 -R200, with a logarithmi c slope 
of ~ 2 - 2.3 in the range (0.3 - 1) feno. lEckert et al.l (120111 ) 
performed a stacking of the density profiles of a sample of 
clusters observed through ROSAT to analyze the outskirts 
of clusters, although they cannot determine any spectral in- 
formation (i.e. the gas temperature) from these data. Their 
average density profile steepens beyond -R500, in agreement 
with the present work and with previous w orks by ROSAT 
|Vikhlinin et al.lll999l ). lEckert et ail (120111 ) also argued that 
the shallow density profiles observed in some clusters by 
Suzaku might be induced by observations in preferential di- 
rections (e.g. along filaments) and do not reflect the typical 
behavior of cluster outer regions. 

It is also interesting to observe that the nor malization of 
the av erage density near the virial radius from lEckert et all 
|201ll ) is about a 50% higher than the present analysis, 
where we found a v ery good agreement wit h the predictions 
from simulations of lRoncarelli et all (|2006l ). Unlike SZ data, 
X-ray brightness is indeed biased by dense, cold clumps in 
the outer volumes, being traced by the square of the gas 
density, and therefore boosting the same gas density. 

We also observe a good agreement of our slope 
of the tempera t ure w ith the theoretical predictions of 
iRoncarelli et all l|2006tl out to 0.7i?2oo, though near -R200 
our temperature profile is steeper, suggesting the presence 
of cold, clumpy gas. 

Now we focus on the entropy profile S = kT/n^ 3 , since 
this is a powerful tool to trace the thermal history of the 
ICM. Indeed the gas entropy records the thermodynamic 
history of the ICM as the product of both gravitational and 
non-gravitat ional proce sses, shaping its observed structure 
accordingly (I Void looHal ) . The measurements of the gas en- 
tropy in the inner regions (0.1 -R200) show that the observed 
value of S is higher than th e expected from the adiabatic sce- 
nario (|Ponman et al.lll999l ) , which include only gravity. This 
excess in the entropy (labeled as entropy "floor" or "ramp" ) 
with respect to the prediction of the adiabatic model calls 
for some energetic mechan ism, in addition to g r avity, such as 
(pre)-heating and coolin g |Borgani et a l. 2005; B ryan ||2000| ; 
iMorandi fe Ettoril l2007l l . Somehow these non-gravitational 
processes intervene to break the expected self-similarity of 
the IC gas in the innermost regions. Nevertheless, in the 



outer volumes simple theoretical models predict that the 
entropy S should be self-similar and behave as a power-law 
with radius. Models of shock dominated spherical collapse 
show that matter is shock heated as it fall s into clusters un- 
der th e pull of gravity, with a slope of ~ 1.1 (|Tozzi fe Norman] 
l200ll ). In the present work, we aim at comparing theoretical 
predictions with the observed entropy profile out to i?200- 

It is interesting to point out that our entropy pro- 
file (see Figure [4]) roughly follows this expected trend for 
0.2 -R200 < R < O.8-R200 with a logarithmic slope of ~ 1. 
In particular, in this spatial range w e find a good ag ree- 
ment with the adiabatic predictions of IVoit et alj (|005bl ) by 
considering SPH simulations through the GAD GET code 
l|Springel et al.ll200ll ; ISpringel fc Hernquistll2 002 ) , where the 
entropy is defined as 



S(r) = S200 1.32(r/r 2 oo) 



(18) 



5*200 being a characteristic value of the entropy at th e over- 
density of 200 (see, e.g., Eq. 2 in IVoit et ail l005bh . Voit 
and collaborators also consider the semi-analytical mod- 
els by using clusters simulated by the AM R code ENZO 
|Norman fc Brvadfl999l ; lO'Shea et aT]|2004l ): in the latter 
case the normalization of the above theoretical relation is 
~ 10 per cent higher than in the SPH simulations, in bet- 
ter agreement with our constraints. This suggests that for 
R > 0.2 -R200 the physics of the X-ray emitting gas are rela- 
tively simple and nearly self-similar, and SZ measurements 
can be used robustly for cosmological works. 

Nevertheless, we observe a flattening of S from the 
power-law shape in the outskirts (R > O.8-R200), as in- 
ferred from Suzaku X-r ay observations 1 George et alj|2009l ; 
ISimionescu et al.ll201 ll ). perhaps indicative of infal ling gas 
which is not dynamically stable (|Nagai fc Laull201ll ). While 
numerical simulations predict gas clumping in the cluster 
outskirts, SZ data are less biased by dense, cold clumps and 
gas in a multi-phase state in the outer volumes with respect 
to X-ray data, since the SZ intensity depends linearly on the 
density, unlike the X-ray flux. Indeed the previous Suzaku 
X-ray measurements point to a more pronounced flattening 
of S in outer regions with respect to the present analysis. 

The gentle flattening of the entropy profile in the out- 
skirts suggests a need for cool phase gas and non-thermal 
pressure support in order t o maintain dy n amic stability, 
as already deduced in i]6.3l lAfshordi et all l|2007l) stacked 
WMAP observations of a large sample of massive clusters 
and found a deficit in thermal energy in the outskirts from 
the Sunyaev-Zeldovich profile, also arguing for a cool phase 
of the ICM. Yet, it is not well understood how these different 
phases would mix, and complicated gas physics would likely 
result, as suggested from the marginal disagreement in the 
level of non-thermal pressure between simulations and the 
current work. 



7 SUMMARY AND CONCLUSIONS 

In this paper we employed a physical cluster model for 
Abell 1835 with a triaxial mass distribution including sup- 
port from non-thermal pressure, proving that it is consis- 
tent with all the X-ray, SZ and lensing observations and 
the predictions of CDM models. This model relies on the 
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following assumptions: 1) the use of a dual Pseudo Isother- 
mal Elliptical Mass Distribution (dPIE) mass model for the 
lensing data; 2) the assumption of generalized hydrostatic 
equilibrium (HE), which also accounts for the non-thermal 
energy contribution; 3) the non-thermal contribution traces 
the thermal pressure, up to a scale factor taken to be a power 
law with the cluster radius; 4) the DM halo density follows 
a generalized Navarro, Frenk and White (gNFW) triaxial 
model. We stress that, given the complementary data sets 
that have been included in this work, we do not need to rely 
on any prior from numerical simulations. 

We presented the first observational measurement of the 
non-thermal pressure out to J?200- The level of non-thermal 
pressure has been evaluated to be a few percent of the to- 
tal energy budget in the internal regions, while it reaches 
about 20% in the outer volumes, a value which is lower than 
the predictions from numerical simulations. This indicates 
that an improved physical treatment in the ICM might be 
needed in simulations to match the amount of non-thermal 
energy with our observational findings. This has important 
consequences for estimating the amount of energy injected 
into clusters from mergers, accretion of material or feedback 
from AGN. 

We analyzed the physical properties of the IC gas in the 
outer volumes out to -R200, focusing on the entropy, which 
is a powerful tool to trace the thermal history of the IC gas. 
We find a good agreement with the theoretical predictions, 
indicating that outside the innermost regions the physics of 
the X-ray emitting gas is relatively simple and nearly self- 
similar. Nevertheless, we observe that entropy tends to gen- 
tly flatten in the outer volumes, likely indicative of infalling 
clumpy and cold gas which is not dynamically stable. In this 
perspective SZ data are crucial for unbiased measurements 
of the cluster physical properties out to large radii, like those 
presented here. 
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